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This note presents an overview of general-purpose and specific Monte-Carlo event generators used in the 
simulation of the CERN - CHORUS experiment, aiming to search for — > v T oscillations and charm particle 
decays in an emulsion target. 



1. Introduction 

CHORUS is an experiment primarily designed 
to search for — » v T oscillations in the CERN 
wide-band neutrino beam. It also aims to investi- 
gate a wide range of charm physics. In four years 
of exposure, more than 10 6 neutrino interactions 
have been accumulated in the emulsion target, 
out of which several thousands lead to the pro- 
duction of a charmed particle. Using its massive 
calorimeter as a target, CHORUS also has a po- 
tential to study other neutrino physics for which 
large statistics is the crucial point. 

Monte-Carlo (MC) simulation is an essential 
part of the experiment. It includes a detailed 
model of the neutrino beam, generation of neu- 
trino events, full GEANT-based simulation of the 
detector response and the reconstruction. 

This report is organized as follows. In Section 

2, a brief description of the CHORUS physics and 
experimental set-up is given. The MC procedure 
is reviewed in Section 3. Section 4 is devoted to 
general-purpose event generators and their verifi- 
cation on the existing experimental data. In Sec- 
tion 5, more specific generators are considered. 



2. CHORUS: physics and set-up 

The CHORUS detector § is a hybrid setup. 
A 770-kg nuclear emulsion is used as the primary 
target for neutrino interactions, allowing three- 
dimensional reconstruction of short-lived parti- 
cles like the r lepton and any charmed hadron. 
The emulsion target consists of four stacks. Each 
stack is followed by three interface emulsion 
sheets and by a set of scintillating fibre tracker 



planes. They provide accurate predictions of par- 
ticle trajectories into the emulsion stack for the 
location of vertex positions. The emulsion scan- 
ning is performed by fully automatic microscopes 
equipped with CCD cameras and a read-out sys- 
tem. The electronic detectors downstream of the 
emulsion target include a hadron spectrometer 
which measures the bending of charged particles 
in an air-core magnet, a lOOt calorimeter where 
the energy and direction of showers are mea- 
sured and a muon spectrometer which determines 
charge and momentum of muons. In addition, the 
calorimeter was used as an active (lead) target for 
neutrino events. 

The West Area Neutrino Facility at CERN pro- 
vides a beam of 27 GeV average energy consisting 
mainly of with a 5% contamination. During 
the four years of operation, emulsion target mate- 
rial was exposed to the beam with an integrated 
intensity which corresponds to 5.06 x 10 19 protons 
on target. The data from the electronic detectors 
were analysed and the events possibly originating 
from the emulsion stacks were identified. 

Since the CHORUS experiment was designed 



primarily to search for v„ 



oscillation, the 



data selection for the first phase of the analysis 
was optimized for the detection of a r decaying 
into a single charged particle. It allowed to reach 
the limit on the — > v T oscillation probability as 
Ppr < 3.4 x 1CT 4 at large Am 2 value @. CHO- 
RUS has also seen neutrino induced diffractive 
D* production || and charm pair production in 
charged-current (CC) interactions In addi- 
tion, with the calorimeter target, neutral-current 
(NC) production of 3 ftp was observed || and pre- 
liminary results on differential cross-sections and 
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structure functions on lead were obtained |6| . 

Recently, the second phase of CHORUS scan- 
ning and analysis has been started. It is based 
on new reconstruction software and improved MC 
packages. The key point is the new emulsion scan- 
ning method called 'netscan'. This method, orig- 
inally developed for the DONUT experiment |t]] , 
consists of recording all track segments within a 
large angular acceptance in a volume surround- 
ing the located vertex position. It can improve 
the limit on the i/„ — > v T oscillation probability 
by a factor of 3 to reach the design goal. In the 
second phase, a search for charm decays in the 
emulsion is also performed. Several thousands of 
such events are expected to be located. Based 
on a small part of the statistics, the D° produc- 
tion cross section has already been measured with 
an accuracy better than 10% ||. The following 
charm physics is of interest: 

• meson rates, branching ratios for different 
channels, extraction of CKM matrix ele- 
ments \V c d\ and \ V CS \, strange sea of the nu- 
cleon and fragmentation functions; 

• charmed baryon production and decays; 

• diffractive production of D* to measure dif- 
ferent branching ratios and determine \ V CS \; 

• associated charm production; 

• rare charm decays. 

With the calorimeter as a target, it is possible to 
measure d 2 a/dxdy and structure functions from 
Vp and interactions on lead, to study dimuon 
production and to search for rare processes. Po- 
tentially the same tasks can be also fulfilled using 
the first two magnets of the iron muon spectrom- 
eter as an active target. In addition, with a set 
of lead, iron, marble and plastic targets[] placed 
in the neutrino beam ratios of CC total cross- 
sections are being determined. 

3. CHORUS MC simulation procedure 

The procedure of MC simulation in the CHO- 
RUS experiment consists of modelling the neu- 

1 At the end of experiment the emulsion target was re- 
placed by this set of 4 targets. 



trino beam, the generation of physical events, the 
simulation of the detector response and the recon- 
struction. 

3.1. Simulation of the neutrino beam 

Parent mesons for neutrinos, modelled with 
FLUKA § by NOMAD are fed into the 

CHORUS GEANT 3.21-based neutrino beam 
simulator called GBEAM O]. It performs track- 
ing of these mesons through the neutrino tunnel 
and models their decays into neutrinos. The out- 
put of the simulation is the flavour, the creation- 
vertex and the four-momentum of the neutrino. 

3.2. Generation of physical events 

Deep-inelastic (DIS) events are simulated with 
the JETTA @ package. (Quasi)-elastic (QE) 
processes and resonance production are modelled 
with the RESQUE Q code. Information about 
these general-purpose packages can be found in 
Section 4. For specific processes, dedicated gener- 
ators were developed (Section 5). All event gen- 
erators use incoming neutrino information pro- 
duced by GBEAM. Their output is similar to 
the LUND LUJETS common block Jl| (particle 
codes, momenta, vertices and history). 

3.3. Simulation of the detector response 

The package which simulates the detec- 
tor response, called Eficass (Emulsion Fibers 
CAlorimeter and Spectrometer Simulation), con- 
tains a very detailed description of the detector 
geometry in the GEANT 3.21 framework. The 
longitudinal vertex positioning is performed with 
a special geantino-based algorithm in accordance 
with the material density the neutrino traverses. 
Eficass digitizes GEANT hits in each active vol- 
ume of the CHORUS detector subsystem. Its 
output, similar to what is produced by the data 
acquisition system, is sent to the reconstruction 
program. 

4. General-purpose generators 

4.1. JETTA 

JETTA (JETs in Tau Analysis) (if] was writ- 
ten for the CHORUS kinematical pre-selection of 
events by embedding into a general DIS genera- 
tor (LEPTO 6.1 @) the formalism of Ref. @ 
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to describe interactions with a massive charged 
lepton (5 structure functions). The code takes 
care of the polarization effects of the outgoing lep- 
ton also with respect to its decays by using the 
TAUOLA Hi package. The rest, is plain JET- 
SET H code, tuned on the BEBC data |||| 
and it deals reasonably well with some charm 
production aspects and different interactions (NC 
and CC). JETTA generates DIS events (W 2 >2 
GeV 2 ) initiated by v e , or v T with energies 
in the region of 3-^300 GeV. Structure functions 
GRV94LO [glj] are used by default, other choices 
are possible via steering cards. Nucleon composi- 
tion of the target and Fermi motion are also taken 
into account. JETTA is able to simulate events of 
specific classes like open charm and dimuon pro- 
duction. Threshold behaviour of the cross section 
due to the non-zero charm quark mass (m c ) is in- 
troduced but without slow rescaling. The Peter- 
son longitudinal fragmentation function ]22] ] with 
e = 0.072 |2^] is used for charm quark hadroniza- 
tion. 

JETTA reproduces reasonably well neutrino 
bubble chamber experimental data on charged 
track multiplicities, transverse shower develop- 
ment, pion energies and fragmentation functions 
for hadrons. ft was also checked on our CC, 
charm and dimuon data obtained both with emul- 
sion and calorimeter targets. 

4.2. RESQUE 

RESQUE @, an event generator for RESo- 
nance and (QUasi)Elastic events, has been de- 
rived from a preexisting code. Computation of 
cross sections, as well as the routines which de- 
scribe resonance decays have been adapted from 
the Soudan-II RSQ generator p4} ]. The main 
changes concern: 

- significant extension of the energy range for the 
calculation of cross-sections and for the genera- 
tion of events; 

- the incorporation of nuclear effects, Fermi mo- 
tion pEf and Pauli suppression; 

- the inclusion of v r interactions, and r decays, 
taking into account the polarization of the lepton 
in the final statej^j 



Three nucleon formfactors are used in the calcu- 
lation of QE CC cross section Pq| . The fourth 
formfactor which is proportional to mass squared 
of the outgoing lepton, is taken into account in 
the case of v T interactions. For elastic NC inter- 
actions, the cross sections are taken from Ref. |27| . 

The original routines from RSQ code are 
used to generate the contribution of 16 baryonic 
resonances (N* and A) with invariant-mass W < 
2 GeV. All resonances and their decays are gen- 
erated separately. The production cross-sections 
have been computed in Ref.|^] using the FKR 
pS| ] semirelativistic quark model. Total cross sec- 
tions of QE vn interactions and A ++ production 
in vp- interactions are found to be in an agree- 
ment with experimental data |50] ^2| . 

4.3. Verification of CHORUS MC chain 

The combined performance of JETTA and 
RESQUE generators was checked on CHORUS 
data. Fig.l shows the reconstructed neutrino 
and antineutrino energy of CC events in the 
calorimeter. The difference between detector- 
deconvoluted data and MC does not exceed 
10-^15%. It is partially due to imperfections in 
our simulation of neutrino parent meson trans- 
portation through the decay tunnel. Charged 
multiplicity of tracks predicted in the emulsion is 
presented in Fig. 2. The agreement between data 
and MCQ is quite satisfactory. 

4.4. MICKEY 

MICKEY §] is used as a fast MC genera- 
tor for CC interactions in the calorimeter. In- 
stead of doing a full GEANT simulation, it ap- 
plies parameterized detector resolution and inef- 
ficiencies. These parameterizations have been ob- 
tained from test-beam data and from full detector 
simulations. The parton distribution set used in 
MICKEY is that of GRV94LO Q. Violation of 
the Callan-Gross relation is modelled using a pa- 
rameterization of R(x, Q 2 ) taken from a fit to the 
world data [[|| . Target mass effects and non-zero 
m c (slow rescaling) are taken into account. The 
structure functions are multiplied by a nuclear 
correction as a function of x according to Ref. Q . 



2 It is assumed that the generated r has always he linit y -1 
as obtained from an extrapolation from DIS case jl7|. 



3 Both data and MC were passed through the same recon- 
struction program. 



4 



QED radiative corrections are calculated follow- 
ing the scheme of Ref.(3f|. The generator was 
tuned and validated using our calorimeter and 
emulsion data f||36|]. 

5. Specific generators 

Four generators were developed for charm stud- 
ies: 

- ASTRA for diffractive D s production jj7); 

- QEGEN for QE production of charmed baryons 

f§; 

- MCDIMUON for dimuon - involving charm - 
production HQ; 

- CCBAR for pair production of charmed parti- 
cles both in CC and NC f§. 

In addition, there is a generator for detailed sim- 
ulation of "white star kinks", WSK Q. 

5.1. ASTRA 

The core of ASTRA Q consists of a set of 
functions F n describing the matrix elements for 
diffractive production of (vector) D* and direct 
production of (pseudoscalar) D s , in neutrino in- 
teractions, according to the choice of theoretical 
models indexed by n. It is left to the user to 
select her/his favourite model. By default, mod- 
els from Ref.Q (Ref.@) are used for D* (D s ) 
production. The first approach is based on gener- 
alized vector dominance while the second one uses 
PCAC as a fundamental ingredient. Non-zero m c 
is taken into account in both. 

It should be noted that a correctly selected in- 
tegration algorithm and correctly computed in- 
tegration limits are crucial. Ref.jli} contains a 
thorough study of physical and computational 
pitfalls to be avoided. 

Fragmentation and decays are performed with 
JETSET H). Nuclear effects other than an ad- 
hoc diffractive slope are not included in ASTRA. 

5.2. QEGEN 

In the QE charm production a valence <i-quark 
is transformed, through weak interaction, into a 
c-quark. The other two valence quarks in the nu- 
cleon behave as spectators, so only one hadronic 
state (charmed baryon A+, £+, £++ or 

is produced and finally each event will 
contain a A+ baryon. There are two classes of 



theoretical models which try to describe the QE 
charm production. The first class is based on 
SU(4) flavour symmetry which is badly broken. 
A different approach is based on local duality 
in i/N scattering on the basis of QCD. Although 
the predictions of the models for the total cross 
sections differ by one order of magnitude, Q 2 - 
distributions are very similar. To simulate the 
dynamics of QE charmed baryon production, the 
differential cross sections from Ref.|4j| are used. 
JETSET is employed to perform fragmenta- 
tion and charmed baryon decays. The model to 
generate Fermi motion in emulsion target nuclei 
is adopted from Ref. |25| . Note that experimental 
data on QE charm production are very poor. 

5.3. MCDIMUON 

The event generator MCDIMUON @ was 
originally developed to simulate opposite-sign 
dimuon events in the CHARM-II detector p3( ]. 
CTEQ3L Q| parton distribution functions are 
used. The masses of the charm quark and the 
nucleon are not neglected in the hard-scattering 
cross-section (slow rescaling is included). The 
fragmentation of charm quarks is described by 
the fragmentation function from Ref. |2^] . Excited 
charmed hadrons are not taken into considera- 
tion. The transverse momentum distribution is 
chosen as an exponential function of p 2 . A cor- 
rection for radiative effects is calculated for the 
inclusive CC cross-section using the prescription 
of Ref. [3f| . Decays of charmed hadrons only de- 
pend on the phase space. The nucleon content of 
the target is also taken into account. 

The code is very flexible and the user can pro- 
vide many parameters through steering cards. It 
was tuned for simulation of dimuon events in the 
CHORUS calorimeter Eol. MCDIMUON was 



verified both on CHARM-II and CHORUS data. 

5.4. CCBAR 

The simulation of charm pair production in 
NC and CC interactions is performed using the 
general-purpose HERWIG event generator |4^| , 
version 6.1 |4^| . Several features of this gener- 
ator, like heavy flavour hadron production and 
decay via QCD coherence effects and the clus- 
ter hadronization of the jets via non-perturbative 
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gluon splitting, are especially relevant. The 
parton-shower approach is used both for ini- 
tial and final states. MRS 5 (Owens) structure 
functions are used jf9|. The Peterson's model 
for fragmentation functions is used with e = 
0.072. Nuclear effects such as evaporation and 
re-interactions are not simulated. 

The process in NC is based on the Z°-gluon 
fusion. A dedicated mode to produce this pro- 
cess is foreseen by the generator. Associated 
charm production in CC is based on the split- 
ting of cc quarks from a gluon which is radiated 
through bremsstrahlung of a light quark. No ded- 
icated mode is foreseen for such a process which 
is therefore obtained among conventional CC in- 
teractions. It should be noted however that total 
cross sections predicted by HERWIG Q and cal- 
culated in Ref. |51j differ by one order of magni- 
tude. The CCBAR generator cannot be validated 
as the available experimental data for associated 
charm production are very poor |^|,[52|,^3| . 

5.5. WSK 

A hadron-nucleus interaction which produces 
only one mimimum ionizing particle and no heav- 
ily ionizing tracks or other signs of nuclear break 
up is usually called "white star kink" (WSK). 
Note that, this is the main background for os- 
cillation searches in CHORUS in "0/x" mode §, 
which is very difficult to estimate. 

The simulation ^lj is composed of two sepa- 
rate parts. The first is the generation of hadron 
interactions in emulsion, the second is the simu- 
lation of emulsion response to the different parti- 
cles. The hadron-emulsion interactions have been 
simulated using FLUKA The emulsion re- 
sponse description is a very delicate point. It is 
based on the range in emulsion, for which the pa- 
rameterization is obtained by an old experiment 
|54| . The generator has been tested and it repro- 
duces the CHORUS data § quite well. 

6. Conclusions 

The CHORUS experiment has two general- 
purpose MC event generators which were vali- 
dated using available experimental data. Tun- 
ing of dedicated generators, developed for charm 



physics, is well advanced. Currently, a complete 
set of CHORUS MC generators is available. 

REFERENCES 

1. CHORUS Collaboration, E. Eskut et al, Nucl. 
Instr. Meth. A401 (1997) 7. 

2. CHORUS Collaboration, E. Eskut et al., 
Phys. Lett. B497 (2001) 8. 

3. CHORUS Collaboration, P. Annis et al. 
Phys. Lett. B435 (1998) 458. 

4. G. De Lellis, Ph.D. Thesis, Universita degli 
Studi di Napoli "Federico II" , Italy, 2000. 

5. CHORUS Collaboration, E. Eskut et al., 
Phys. Lett. B503 (2001) 1. 

6. R.G.C. Oldeman, PhD thesis, University of 
Amsterdam, The Netherlands, 2000. 

7. DONUT Collaboration, K. Kodama et al., 
Phys. Lett. B504 (2001) 218. 

8. CHORUS Collaboration, A. Kayis-Topaksu 
et al., preprint CERN-EP/2002-005, accepted 
for publication in Phys. Lett. B. 

9. A. Fasso et al., SARE-3 Workshop, KEK Re- 
port Proceedings 97-5, p.32 (1997). 

10. A. Guglielmi, private communications. 

11. GEANT 3.21, CERN program library long 
write up W5013. 

12. S. Sorrentino, Diploma Thesis, Naples Uni- 
versity, Italy (1995). 

13. P. Zucchelli, PhD thesis, Universita di Fer- 
rara, Italy (1995). 

14. S. Ricciardi, Ph.D. thesis, Universita Ferrara, 
Italy (1996). 

15. T. Sjostrand, Comp. Phys. Comm. 82 (1994) 
74. 

16. G. Ingelman, preprint TSL/ISV 92-0065, Up- 
psala University, Sweden (1992). 

17. C.H. Albright and C. Jarlskog, Nucl. Phys. 
B84 (1975) 467. 

18. S. Jadach et al., Comput. Phys. Commun. 76 
(1993) 361. 

19. BEBC Collaboration, H. Grassier et al., Nucl. 
Phys. B214 (1983) 369, Nucl. Phys. B223 
(1983) 269. 

20. D. Allasia et al., Z. Phys. C24 (1984) 119. 

21. M. Gliick, E. Reya and A. Vogt, Z. Phys. C67 
(1995) 433. 

22. C. Peterson et al., Phys. Rev. D27 (1983) 105. 



G 



23. 

24. 

25. 

2G. 

27. 

28. 

29. 

30. 
31. 
32. 
33. 

34. 
35. 

3G. 

37. 



38. 
39. 

40. 

41. 

42. 

43. 

44. 

45. 

4G. 

47. 

48. 
49. 



CHARM-II Collaboration, P. Vilain et al., 
Europ. Phys. J. Cll (1999) 19. 
G.D. Barr, Ph. D. Thesis, Keble College, Ox- 
ford, Great Britain (1987). 
A. Bodek and J.L. Ritchie, Phys. Rev. D23 
(1981) 1070. 

C. H. Llewellyn-Smith, Phys. Rep. 3 (1972) 
261. 

R.E. Hendrick and L.-F. Li, Phys. Rev. D19 
(1979) 779. 

D. Rein and L.M. Seghal, Annals of Physics 
133 (1981) 79. 

R.P. Feynman, M. Kislinger and F. Ranvdal, 

Phys. Rev. D3 (1971) 2706. 

J. Bell et al., Phys. Rev. Lett. 41 (1978) 1008. 

P. Allen et al., Nucl. Phys. B176 (1980) 269. 

S. Belikov et al., Z. Phys. A320 (1985) 625. 

L.W. Whitlow et al., Phys. Lett. B250 (1990) 

193. 

J. Gomez et al., Phys. Rev. D49 (1994) 4348. 
D.Yu. Bardin and V.A. Dokuchaeva, Preprint 
JINR E2-86-260, 1986. 

O. Melzer, Ph.D. thesis, University of Ams- 
terdam, The Netherlands, 2001. 
O. Melzer, Diploma Thesis, Westfahsche 
Wilhelms-Universitat, Minister, Germany, 
1997. 

F. Di Capua, private communications. 

V. Lemaitre, PhD Thesis, University 
Catholique de Louvain, Belgium, 1995. 

C. A.F.J. van der Poel, Ph.D Thesis, Univer- 
sity of Nijmegen, The Netherlands, 1999. 

A. Satta, Ph. D. Thesis, Universita degli 
Studi di Roma "La Sapienza", Italy, 2001. 
M.K. Gaillard, S.A. Jackson and D.V. 
Nanopoulos, Nucl. Phys. B102 (1976) 326. 

D. Rein and L.M. Seghal, Nucl. Phys. B223 
(1983) 29. 

M.S. Chen, F.S. Henyey and G.L. Kane, Nucl. 
Phys. B118 (1977) 345. 

R.E. Shrock and B.W. Lee, Phys. Rev. D16 
(1976) 2539. 

CTEQ Collaboration, H.L. Lai et al., Phys. 
Rev. D51 (1995) 4763. 

G. Marchesini et al., Comp. Phys. Comm. 67 
(1992) 465. 



Phys. Lett. B306 (1993) 145, erratum Phys. 
Lett. B309 (1993) 492. 

50. M.H. Seymour, Nucl. Phys. B436 (1995) 163. 

51. K. Hagiwara, Nucl. Phys. B173 (1980) 487. 

52. NuTeV Collaboration, A. Alton et al., Phys. 
Rev. D64 (2001) 012002. 

53. E531 Collaboration, N. Ushida et al., Phys. 
Lett. B206 (1988) 375. 

54. W.H. Barkas, "Nuclear research emulsions", 
Academic Press, 1963. 



^0000 

Q5000 
ZJ 

£20000 










































o 

^3 5000 






















'J5 






















^0000 
5000 



























4 


6 


8 


1 


1 


1 


1 


1 


2C 



Energy of neutrino, GeV 




Energy of ontineutrino, GeV 



Figure 1. Neutrino (top plot) and antineutrino 
(bottom plot) energy spectra for CC events in 
the calorimeter. The solid histograms are for ex- 
perimental data, the dashed histograms for MC. 
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Charged multiplicity 

Figure 2. Charged multiplicity for events having 
predictions of particle trajectories into the emul- 
sion stack for the location of the vertex positions. 
The solid histogram is for experimental data, the 
dashed histogram for MC. 
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